Re-assigning (1x2) reconstruction of rutile TiO2(110) from DFT+ U calculations 
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Physically reasonable electronic structures of reconstructed rutile TiO2(110)-(lx2) surfaces were 
studied using density functional theory (DFT) supplemented with Hubbard U on-site Coulomb 
repulsion acting on the d electrons, so called as the DFT+[7 approach. Two leading reconstruction 
models proposed by Onishi-Iwasawa and Park et al. were compared in terms of their thermodynamic 
stabilities. 
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Rutile Ti02 and its surfaces represent model systems 
to explore the properties of transition metal oxides that 
are important in technological applications such as catal- 
ysis, photovoltaics, and gas sensing [l|, to name a few. 
Truncated or stoichiometric (110) surface of rutile is the 
most stable one among all surfaces of titania Q. Upon 
thermal annealing or ion bombardment TiO2(H0)-(lxl) 
surface is reduced by loosing the bridging oxygens, and 
is often undergo a (1x2) reconstruction with row for- 
mations SHU. The identification of these rows on re- 
constructed surfaces, in three dimensions, is difficult by 
experimental methods [Hj. The best candidate for mod- 
eling this reconstruction involves the addition of a Ti20s" 
molecule on the surface unit cell (added-row model) pro- 
posed by Onishi and Iwasawa [3] . In addition to theoret- 
ical studies it was supported by electron stimulated des- 
orption of ion angular distribution (ESDIAD), scanning 
tunnelling microscopy (STM), and low-energy electron 
diffraction (LEED) experiments [5l4ll). On the other 
hand, Shibata et al. [13], using advanced transmission 
electron microscopy (TEM) observations reported results 
that are consistent with a different model, proposed by 
Park et al. [lH, for which the additional unit is "T12O". 
The main difference between the Onishi-Iwasawa and 
Park et al. models (Onishi and Park models, respec- 
tively, in what follows) is the locations of Ti interstitial 
sites on the surface [l3[ (see Fig. [I]). 

Now there is a debate about which formation gives rise 
to the (1x2) long range order, the last proposed one or 
the best previous candidate? In this context we study 
the electronic properties of these two models using Hub- 
bard U corrected total energy density functional theory 
(DFT+Z7) calculations to get physically reasonable re- 
sults comparable to existing experimental data. We dis- 
cuss which of these leading models can be assigned to 
describe the (1x2) reconstructed surface by comparing 
them according to their thermodynamic stabilities. 

Band structures of reduced and reconstructed 
TiO2(110) surfaces determined by pure DFT calculations 
does not agree with experiments Io|, [HI, Ell • Failure of 
the standart DFT is not limited to band-gap underesti- 
mation stemming from the many-electron self-interaction 




* UJXG 

» i 1 y 




FIG. 1. Fully relaxed (1x2) reconstructed rutile TiO 2 (H0) 
surface models of (a) Onishi-Iwasawa 3] and (b) Park et 
al. [13] showing 3 Ti- layers from the top. 



error (SIE). More importantly, it does not predict exper- 
imentally observed gap states that are associated 
with the excess electrons due to the formation of surface 
oxygen vacancies. In DFT calculations, these Ti 3d elec- 
trons occupy the bottom of the conduction band (CB) 
giving metallic character. Hence, hybrid DFT meth- 
ods need to be used. For instance, SIE can be partly 
corrected by partially mixing nonlocal Fock exchange 
term with DFT exchange term [H, [HI or DFT+/7 ap- 
proach [20[ can make up for the lack of strong correla- 
tion between the 3d electrons, a shortcoming of common 
exchange-correlation functionals. Empirical Hubbard U 
term accounts for the on-site Coulomb repulsion between 
the Ti 3d electrons. By examining different values of the 
U parameter, experimentally observed gap state of the 
reduced TiO gCllO) surface was obtained ~0.7-0.9 eV be- 
low the CB [2lM23|. 

For the (1x2) reconstructed surface with Ti203 added 
row, Kimura et al. Q, and then, Blanco- Rey et al. [To[, 
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obtained the Ti 3d states positioned inside the CB, and 
proposed this model to be metallic by pure DFT meth- 
ods. Recently, using spin polarized DFT+Z7 calculations 
with a suitable choice of the U parameter, we have shown 
that the Onishi model can be semiconducting [25j a s 
the reconstructed surface is observed experimentally [26| . 
Important questions still remain to be answered such as 
where the excess charge density is distributed and which 
model structure describes the (1x2) reconstruction. 

We used Perdew-Burke-Ernzerhof (PBE) [27} gra- 
dient corrected exchange-correlation functional supple- 
mented with Dudarev's U term [io[ as implemented in 
the VASP code [Hj . The ionic cores and valence electrons 
were treated by the projector-augmented waves (PAW) 
method 29|, 30[ up to a cutoff value of 400 eV. In order 
to get well converged energetics, we adopted stoichiomet- 
ric slabs with 10 Ti-layers (30 atomic layers) separated 
by -15 A vacuum from their periodic images. We built 
the Onishi and the Park models of (1x2) reconstructed 
rutile (110) surfaces by adding Ti203 and T^O groups, 
respectively, along [001] both at the top and at the bot- 
tom of the slabs (see Fig. [I]). Hence, having an unphysical 
dipole across the slab and getting two different groups of 
surface states in the gap has been avoided. Instead, sym- 
metric slabs bring about the same group of surface states 
(degenerate) from both surfaces. 

Our choice of U=5 eV follows from our examination 
of the effect of different U values on both the geometry 
and the energy bands of rutile Ti02- While a larger U 
value can give a wider band gap, it would make a sig- 
nificant distortion in the atomic structure. Inclusion of 
Dudarev U =5 eV term acting on Ti 3d electrons gave rea- 
sonable values for the atomic positions with bond lengths 
in agreement with the experimental data for the stoichio- 
metric, reduced, and the reconstructed rutile (110) sur- 
faces. This choice is also consistent with the previous 
calculations describing the reduced 21-23| and recon- 



TABLE I. PBE+[7 results of atomic positions for the Onishi 
model with U =5 eV. 



structed [25| surfaces. Spin polarization was also found to 
be important in determining the semiconducting ground 
state of reconstructed surfaces and in correctly describ- 
ing Ti defect states while it is insignificantly small for 
the stoichiometric surface. Our PBE+LT calculations de- 
termined the ground states with spin multiplicities of 2 
and 6 per surface of (1x2) unit cells for Onishi and Park 
models, respectively. 

Calculated atomic coordinates and the LEED re- 
sults [Io|, El are only slightly different as shown in Ta- 
ble U and Table El for each of the model structures. The 
most noticeable deviation is seen in the z component of 
0(4) in the Onishi model. While keeping distortion to 
the atomic positions small, PBE+Z7 with U=5 eV repro- 
duces the gap states as shown in Fig. O Defect states 
has been found to be 1.24 eV and 0.68 eV below the CB 
for Onishi and Park models, respectively, in agreement 
with experiments [iBI, - In the Onishi model, the gap 
state with a dispersion width of 0.65 eV reveal charge 
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LEED data in Ref.pjJ. 



TABLE II. PBE+[7 results of atomic positions for the Park 
model with [7=5 eV. 
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LEED data in Ref.pjJ. 



derealization to the interstitial Ti atom just below in- 
plane Ti5c showing d z i character and to the oxygen atom 
beneath it as shown in Fig. [3l 

As finding the almost correct position of the experi- 
mentally observed Ti 3d states, describing the origin of 
those states is also important. According to the ultra- 
violet photoelectron spectroscopy (UPS) and LEED ex- 
periments, the excess electrons upon O loss are believed 
to delocalize around the surface Ti atoms, adjacent to 
the vacancies [17]. Although this vacancy model is sup- 
ported by some calculations [ljj [2l[, there are studies 
that indicate the role of subsurface Ti atoms about the 
derealization of the excess charge 0, [HI, [HI . Since pure 
DFT gives a clean band gap, based on their STM and 
UPS measurements Wendt et at. proposed that the gap 
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FIG. 2. Density of states (DOS) plots and energy band dia- 
grams for the majority electrons belonging to Onishi [3] and 
Park models of TiO 2 (H0)-(l x2) surface. 



states originate from Ti atoms diffused into interstitial 
sites, not from the surface Ti atoms adjacent to bridging 
O vacancies (32|. Without a need to such an additional 
interstitial Ti atom, derealization of the excess charge to 
the subsurface Ti, responsible for the gap state, emerges 
from DFT+/7 calculations as shown in Fig. [3] 

Surface free energy is a good measure to compare the 
thermodynamic stability of model row formations. Sto- 
ichiometric cell is a stack of Ti408-(lx2) units that act 
as bulk layers around the central regions as taken into 
account in Ref. Our ten Ti-layer slab model can be 
expressed as Ti n 02 n with n = 40. Reduced reconstruc- 
tions are led by surface O removal from the stoichiomet- 
ric surface. For our symmetrical slab cell, Ti n 02 n -m 
represents both the Onishi (m=2) and the Park (m=6) 
models. Then, the surface energies were calculated by 
the relation, 



(7 = 



2A 



( E slab _ nE bu\k + mE Q^ ^ 



where A is the surface unit cell area, £^ lab , £; bulk , and 
E° are the energies of a Ti n 02n-m slab, of a bulk T14O8 
unit (—89.3 eV) and of an oxygen atom in its molecular 
form, respectively. Division by two is because the slab 
has two reconstructed surfaces on its both faces. In an 



experiment, the surface layer is an interface between O2 
gas phase and Ti02 bulk crystal. Thermal equilibrium 
can exist if the chemical potential of the atomic species 
are equal in all these phases that come into contact with 
each other. Therefore, determination of the chemical po- 
tential of oxygen atom, E° , limits the accuracy of the 
calculated surface energies. We found the binding energy 
of an O2 molecule to be 6.07 eV which is significantly 
larger than the experimental value of 5.26 eV [33]. The 
tendency of DFT to overestimate it, was also reported 
previously 21, 34|. Therefore, in order to obtain a more 
reasonable energy value for E° , we adopted the experi- 
mental binding energy of O2 (5.26 eV) and DFT result of 
an isolated O atom (-1.89 eV). It gives -4.52 eV for E° . 
By using this reference chemical potential we calculated 
the formation energy of bulk Ti02 from metallic bulk 
Ti and O2 molecule as 9.72 eV, in excellent agreement 
with the experimental value of 9.73 eV [33|. Therefore 
thermodynamic equilibrium between the surface and the 
bulk crystal can also be reached. 
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FIG. 3. Partial charge densities of the gap states. 

In Fig. E] we compare relative stabilities of stoichiomet- 
ric and reconstructed slabs through their surface energies 
with varying number of Ti-layers. Surface energetics of 
the stoichiometric slab at the Hartree-Fock and at the 
standard DFT levels were reported to exhibit odd-even 
oscillations with the number of layers 33-l37| . Our results 
show that the oscillation of both the stoichiometric and 
reduced surface energies with slab thickness settles down 
by the inclusion of U=b (Fig. [4j). Using ten Ti-layer cells 
we calculated well converged surface free energies as 0.74, 
2.03, and 4.82 J-m -2 , for the stoichiometric and recon- 
structed (Onishi and Park models) surfaces, respectively. 
Our results are significantly different from the calcula- 
tions carried out with pure functionals. For instance, 
Morgan et al calculated the surface energy of stoichio- 
metric case to be 0.58 J-m -2 using GGA and to be 0.83 
J-m~ 2 using GGA+£7 ([7=4.2 eV) over (4x2) supercell 
with five Ti-layers [2l[. The latter is in good agreement 
with our GGA+U value of 0.86 J-m -2 calculated with 
5-Ti-layer stoichiometric (1x2) supercell. For the Ti203 
added row model, Elliot et al. found the surface energy 
using spin-polarized DFT to be 3.29±0.08 J-m -2 which is 
corrected by the absolute energy value of an isolated oxy- 
gen atom [7[. On the other hand, for recent Ti20 model 



of Park et al, no calculations for the surface energy were 
reported. 
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FIG. 4. Calculated free energies of TiO2(H0) surface models. 

Relaxed bulk termination appears to be the most sta- 
ble surface. Once the TiO2(110) surface is reduced (fol- 
lowed by a reconstruction), O2 exposure can not restore 
the stoichiometric form again [32|. Thus, comparison 
between the reconstruction models is more meaningful. 
Formation energy of the Ti203 added row proposed by 
Onishi and Iwasawa (with a ground state spin polariza- 
tion of fi = 2) is only 1.29 J-m -2 higher relative to that 
of the stoichiometric surface. 

Standard DFT incorrectly predicts Ti 3d excess charge 
to occupy the bottom of the CB leading to metalliza- 
tion for the reduced and reconstructed surfaces. PBE+£7 
method with U=5 reproduces experimentally observed 
gap states for the (1x2) reconstructions as well as for 
the oxygen vacancies on the rutile (110) surfaces. Fi- 
nally, inclusion of a suitably chosen U parameter in the 
calculations for the reconstructed rutile TiO2(110) sur- 
face, is a simple and promising way of restoring its semi- 
conducting nature by reproducing the band-gap states 
which arise from derealization of Ti 3d excess charge to 
subsurface Ti sites. According to their surface energies, 
Onishi's added row model is more stable than the Park 
model. Therefore, Ti2C>3 added row model confirms ex- 
isting experimental observations and can still be assigned 
as the (1x2) long range order on the rutile (110) surface. 
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